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INTRODUCTION 

Many  fishes  that  swim  with  paired  pectoral  fins  use  fin  stroke  parameters  that  produce  thrust 
force  from  lift  in  powering  their  underwater  flight.  Flapping  locomotion  mechanisms  are  of 
interest  to  behavioral  biologists,  biomechanics  researchers,  and  engineers  attempting  to  develop 
systems  that  can  match  the  performance  of  the  living  creatures.  One  question  that  arises  in 
investigating  the  mechanisms  responsible  for  the  high  thrust/lift  generation  in  flapping 
locomotion  is  the  importance  of  chordwise  and  spanwise  deformation.  Flapping  foil  propulsion 
has  also  received  considerable  attention  in  the  past  few  years  as  an  alternative  to  the  propeller. 
This  mode  of  propulsion  which  involves  no  body  undulation,  has  many  applications,  such  as 
submersibles  propulsion,  maneuvering,  position-keeping  and  flow  control  which  are  of  interest 
to  the  Unmanned  Underwater  Vehicle  (UUV)  hydrodynamic  community  and  unconventional 
aerodynamics  of  Micro  Aerial  Vehicles  (MAV)  for  the  aerodynamic  community.  For  both  types 
of  vehicles,  it  is  important  to  accurately  quantify  the  performance  benefits  of  flexible  wings, 
since  providing  the  flexibility  is  a  major  technical  challenge. 

Dickinson  et  al.  (1999)  has  studied  the  effects  of  the  wing  rotation  in  the  fruitfly,  Drosophila, 
and  Walker  and  Westneat  (1997)  have  studied  the  kinematics  of  the  fin  motion  in  a  class  of 
fishes,  namely  the  bird  wrasse,  experimentally.  Three  dimensional  unsteady  flow'  computations 
over  insects  have  been  carried  out  by  Liu  and  Kawachi  (1998)  who  obtained  qualitative 
agreement  for  the  flow  patterns  over  a  hovering  hawkmoth  with  the  windtunnel  visualizations  of 
Willmott  and  Ellington  (1997).  Ramamurti  and  Sandberg  have  studied  the  force  production  in 
Drosophila  (2002a)  and  in  a  swimming  bird  wrasse  fish  (2002b)  computationally  and  obtained 
good  agreement  with  experimental  results  with  measured  force  data  of  Dickinson  et  al.  (1999) 
and  the  fish  acceleration  data  of  Walker  and  Westneat  (1997).  Ramamurti  and  Sandberg  (2004) 
have  investigated  the  fluid  dynamics  underlying  the  generation  of  forces  during  pectoral  fin 
oscillation  as  fin  rigidity  is  varied. 

The  results  of  that  work  indicated  a  reduced  number  of  ribs  could  be  used  without  losing 
flexibility.  The  primary  objectives  of  this  study  are  to  select  the  shape  of  the  fin  that  allows 
flexibility  while  keeping  the  mechanical  system  less  complex  and  to  select  the  parameters  that 
govern  the  kinematics  of  the  motion  of  the  fin  that  maintains  thrust  for  a  vehicle  in  a  3kt  current 
while  maintaining  the  depth  within  a  specified  range. 

THE  INCOMPRESSIBLE  FLOW  SOLVER 

The  governing  equations  employed  are  the  incompressible  Navier-Stokes  equations  in 
Arbitrary  Lagrangian-Eulerian  (ALE)  formulation  which  are  written  as 


—  +  v.  ■  Vv  +  S7p  =  V-ct, 


0) 

(2) 
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where  p  denotes  the  pressure,  vQ  =  v  -  w  the  advective  velocity  vector,  where  v  is  the  flow 
velocity  and  w  is  the  mesh  velocity  w  and  the  material  derivative  is  with  respect  to  the  mesh 
velocity  w.  Both  the  pressure  p  and  the  stress  tensor  a  have  been  normalized  by  the  (constant) 
density  p  and  are  discretized  in  time  using  an  implicit  time  stepping  procedure.  Thus  the 
equations  are  Eulerian  for  zero  mesh  velocity  and  Lagrangian  if  the  mesh  velocity  is  the  same  as 
the  flow  velocity.  The  present  time-accurate  flow  solver  is  discretized  in  space  using  a  Galerkin 
procedure  with  linear  tetrahedral  elements.  The  details  of  the  flow  solver  have  already  been 
discussed  extensively  elsewhere  (Ramamurti  et.  al.  1992,  1994,  1995,  1999)  in  connection  with 
successfully  validated  solutions  for  numerous  2-D  and  3-D,  laminar  and  turbulent,  steady  and 
unsteady  flow  problems. 


RESULTS  AND  DISCUSSION 

Ramamurti  and  Sandberg  (2004)  have  studied  the  effect  of  variable  rigidity  of  the  Bird- 
wrasse  fin  on  the  fluid  dynamics  of  force  production.  The  computations  showed  that  when  the 
fin  is  made  rigid  by  specifying  the  motion  with  just  the  leading  edge  of  the  fin  tip,  the  thrust 
produced  during  the  upstroke  is  less  than  half  of  the  peak  thrust  produced  by  the  flexible  cases. 
During  the  downstroke,  the  rigid  fin  and  the  fin  with  the  motion  prescribed  with  only  the  leading 
and  trailing  edges  produced  no  positive  thrust,  while  all  the  flexible  cases  considered  reproduced 
the  thrust  production  of  the  fully  deformable  fin.  In  the  case  of  the  rigid  fin,  there  is  a  substantial 
penalty  in  lift  during  the  upstroke. 

Design  and  Construction 

Based  on  this  study,  we  have  reduced  the  number  of  ribs  for  the  mechanical  fin  from  14 
spines  in  the  fish  to  5  ribs,  so  as  to  reduce  the  complexity  of  the  system  while  maintaining 
enough  flexibility  of  the  fin  surface.  The  fin  was  attached  to  a  notional  2ft  long  Unmanned 
Underwater  vehicle  (UUV)  having  a  4"  diameter  cylindrical  cross-section  (Fig.l),  at  an  angle  of 
attack  of  ar  at  the  root  of  the  fin.  The  fin  is  flapped  at  a  frequency  /with  an  amplitude  <|)mlu  about 
the  root  chord. 


Fig.  1.  Schematic  of  the  notional  UUV  with  flapping  fin. 
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The  root  section  of  the  fin  was  selected  to  be  of  rectangular  cross  section  with  rounded 
leading  edges  and  a  tapered  trailing  edge  in  order  to  accommodate  the  actuators  at  the  base  of  the 
ribs.  The  lengths  of  the  ribs  were  obtained  from  the  pressure  distribution  on  the  fin  and  from  the 
structural  analysis  performed  by  UM  researchers.  The  resulting  fin  structure  is  shown  in  Fig. 2. 
The  deforming  motion  of  the  fin  is  obtained  by  prescribing  the  motion  of  the  5  control  points  or 
the  tips  of  the  ribs  (1-5).  The  kinematics  of  the  ribs  is  decomposed  into  two  rotations.  The  first 
one  is  a  rotation  of  the  all  the  ribs  about  the  root  chord  section  termed  as  a  bulk  or  hinge  rotation, 
Fig.  3a,  and  the  second  is  a  supplemental  rib  rotation  which  is  crucial  for  the  deforming  shape  of 
the  fin,  shown  in  Fig.  3b.  For  all  the  results  described,  first  a  steady  state  solution  is  obtained 
using  the  incompressible  flow  solver  described.  From  this  steady  state  as  initial  condition, 
unsteady  flow  over  several  flapping  cycles  is  computed.  The  motion  of  the  fin  is  prescribed 
from  the  kinematics  of  the  5  control  points.  The  resulting  instantaneous  pressure  distribution  is 
integrated  over  the  fin  and  the  body  surface  to  yield  the  thrust  and  lift  forces. 


Fig.  2.  Schematic  of  the  deforming  fin  showing  the  root  cross-section  and  the  5  ribs. 


time  (secs.)  time  (secs.) 

Fig.  3.  Rotation  of  the  hinge  and  the  additional  rotation  of  the  ribs,/=  3.333Hz. 


4 


Ravi  Ramamurli 


Unsteady  Computations 

Initially,  the  fin  was  oriented  to  be  at  an  angle  ar  =  30°  and  the  flapping  amplitude  was  set  to 
76°.  The  position  of  the  beginning  of  the  downstroke  is  obtained  by  appropriately  scaling  the 
amplitude  of  the  flapping  with  that  of  the  Bird-wrasse  fish  and  is  set  to  be  at  -45°  about  the 
root/hinge  axis.  This  configuration  was  first  tested  for  hovering  with  the  fin  flapping  at 
frequencies  of  0.333Hz  and  1.0Hz.  The  frequencies  and  the  amplitudes  were  reduced  from  the 
fish,  which  are  3.333Hz  and  1 14°,  so  that  the  artificial  muscles  will  be  able  to  generate  enough 
force  and  respond  fast  enough.  The  time  history  of  forces  for  these  cases  is  shown  in  Fig.  4  and 
5.  At /=  0.333  Hz,  the  mean  thrust  produced  is  very  nearly  zero  (=-0.008  N),  shown  in  Fig.  4a, 
implying  that  the  fin  produces  nearly  enough  thrust  for  hovering.  Also  a  mean  positive  lift  is 
produced  at  this  frequency.  As  the  frequency  is  increased  to  1  Hz,  a  mean  positive  thrust  of 
0.0 12N  and  the  mean  lift  is  0.05N.  At  hovering,  as  the  frequency  is  increased  the  thrust 
production  during  the  upstroke  increases  sharply.  Although  the  mean  lift  increases  from  0.03N 
to  0.05N,  the  range  of  lift  increases  almost  five  fold  with  increase  in  frequency. 


0.00  3.00  6.00  9.00  12.00  0.00  3.00  6.00  9.00  12.00 

lime  (secs.)  lime  (secs.) 


Fig.  4.  Thrust  and  Lift  production  for  a  hovering  UUV,/=  0.333Hz. 

The  operational  requirements  for  the  UUV  were  set  so  as  to  maintain  position  at  a  3kt  current. 
Computations  were  carried  wit  incoming  velocity  set  to  3kts.  At/=  0.333Hz,  the  current  fin 
configuration  produced  a  drag  of  1.342N  and  a  mean  lift  of  4.023N.  The  flapping  frequency  was 
then  increased  to  1.0  Hz  and  the  mean  drag  and  lift  remained  nearly  unchanged.  The  frequency 
was  further  increased  to  3.333Hz  and  the  mean  drag  reduced  to  1.1 6N  and  the  mean  lift 
increased  to  4.78N.  The  time  history  of  the  force  production  is  shown  in  Fig.  6.  At  this  flow 
speed,  the  current  fin  configuration  and  the  kinematics  is  not  able  to  produce  enough  thrust  to 
maintain  position. 
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Fig.  5.  Thrust  and  Lift  production  for  a  hovering  UUV,/=  1 .0Hz. 


time  (secs.)  time  (secs.) 

Fig.  6.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,/  =  3.333Hz. 

One  of  the  main  differences  between  the  current  fin  design  and  the  wrasse  fin  is  that  the  ribs  in 
the  present  are  parallel  to  each  other.  This  design  was  chosen  to  minimize  the  mechanical 
complexity.  Computationally  the  ribs  can  be  spread  out  without  much  difficulty.  Therefore,  the 
trailing  edge  rib  (rib  5)  was  spread  out  by  60°  with  respect  to  the  first  or  leading  edge  rib. 
Computations  were  carried  out  at  0.333Hz  and  the  mean  drag  remained  nearly  unchanged  and  is 
1.346N,  while  the  mean  lift  increased  to  7.14N.  Due  to  mechanical  complexity  of  spreading  the 
ribs,  this  option  was  put  aside. 

Next,  the  kinematics  of  the  tips  of  the  leading  edge  rib  (ribl)  was  modified  from  the  bird- 
wrasse.  Two  different  methods  were  attempted  to  match  the  kinematics  of  the  fish.  First,  the 
spatial  slope  at  the  leading  edge  fin  tip,  point  1  in  Fig2,  was  matched  to  that  from  the  wrasse. 
The  slopes  at  the  other  control  points,  2-5,  were  retained  from  the  present  kinematics.  This  did 
not  result  in  any  improvement  in  the  thrust.  Next,  the  spatial  slope  at  control  points  2-5  were 
smoothed.  Again,  the  overall  men  thrust  was  still  negative.  The  second  method  is  to  match  the 
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temporal  derivative  of  the  leading  control  point  to  that  of  the  wrasse  with  appropriate  scaling  for 
the  reduced  frequency  of  0.333Hz.  This  also  did  not  result  in  an  improvement  in  the  mean 
thrust.  Few  other  attempts  to  reproduce  some  aspects  of  the  wrasse  kinematics  such  as  the 
maintaining  the  ar  =  53°  or  combining  a  pitching  motion  of  the  fin  in  addition  to  the  flapping  did 
not  result  in  increased  thrust.  These  modifications  can  be  revisited  later. 

Next,  the  ar  was  changed  form  30°  to  0°.  For  this  configuration,  at /  =  0.333Hz,  the  mean 
thrust  improved  from  -1.342N  to  -0.7N,  Fig.  7a,  but  the  mean  lift  suffered  and  reduced  from 
4.023N  to  -1.936N,  Fig.  7b.  As  the  /increased  to  1  Hz,  the  mean  thrust  remained  nearly  the 
same  at  -0.65N  and  with  further  increase  to  3.333Hz,  the  mean  thrust  improved  to  -0.3 IN, 
Fig. 8a,  while  the  mean  lift  was  nearly  -2. ON,  Fig.  8b.  The  loss  of  lift  is  due  to  the  lack  of 
positive  angle  of  attack  of  the  fin. 


Another  key  parameter  that  is  different  from  the  bird-wrasse  is  the  amplitude  of  oscillation, 
which  was  set  to  76°  compared  to  a  valued  114°  in  the  fish.  Hence,  the  amplitude  was  matched 
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to  that  of  the  fish  by  increasing  the  rotations  of  all  the  ribs.  The  combined  rotation  of  the  ribs 
and  the  hinge  for  this  case  is  shown  in  Fig.  9.  Also,  the  angle  of  attack  at  the  root  was  increased 
to  10°.  The  position  of  the  beginning  of  the  downstroke  was  set  to  be  at  -85°  about  the  root- 
chord  axis.  At  a  frequency /=  1.0  Hz,  this  configuration  produced  a  mean  drag  of  0.33N  and  a 
mean  lift  of  -1.476N.  As  the  frequency  is  increased  to  3.333Hz,  a  positive  mean  thrust  of 
0.264N  was  produced,  but  the  mean  lift  remained  negative.  The  results  are  shown  in  Fig.  10. 


time  (secs.) 

Fig.  9.  Rotations  of  the  ribs  for  increased  amplitude, /=  3.333Hz. 

As  the  ar  is  further  increased  to  20°,  the  thrust  and  lift  generation  continued  to  improve.  At/  = 
1.0Hz,  the  mean  thrust  produced  is  -0.25N  and  a  mean  lift  of  0.37N  was  produced,  as  shown  in 
Fig.  11.  At  the  higher  frequency  of  3.333Hz,  a  positive  mean  thrust  and  lift  of  0.33N  and  0.42N, 
respectively,  were  produced,  as  shown  in  Fig.  12a  and  b. 


Fig.  10.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  10°,  (j)max  =  1 14°,/=  3.333Hz. 


As  the  root  angle  of  attack  is  increased  further  to  30°,  the  mean  thrust  reduced  to  -0.7N  at /=  1 
Hz  and  to  -0.2N  at /=  3.333Hz.  The  mean  lift  was  nearly  2.3N  for  both  the  frequencies  tested. 
Therefore,  it  was  deemed  20°  is  near  optimal  angle  of  attack  of  the  root.  In  order  to  produce  a 
prototype  of  the  fin,  it  was  necessary  to  increase  the  spacing  between  the  ribs  from  0.8cm  to  at 
least  1.2cm. 
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Fig.  11.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  20°,  (j)max  =  1 14°,/=  1.0Hz. 

The  spacing  between  the  ribs  was  increased  to  2.0cm.  At  /  =  1.0Hz,  the  mean  thrust  was 
reduced  from  -0.25N  to  -0.63N.  In  an  attempt  to  recover  this  lost  thrust,  the  length  of  the  ribs 
was  doubled.  This  resulted  in  a  positive  thrust  of  0.49N  and  the  mean  lift  of  1.49N,  Fig.  13.  A 
closer  look  at  lift  production  showed  that  the  vehicle  would  undergo  large  excursions  in  the 
vertical  direction.  The  minimum  and  maximum  values  of  the  lift  for  this  case  were  -14.27  N  and 
+  14.04N  compared  to  a  value  of  -0.95N  and  1.91N  for  the  case  when  the  rib  spacing  was  0.8cm 
and  with  the  original  length  of  the  ribs.  Therefore,  the  rib  spacing  was  reduced  to  minimum 
possible  value  of  1.2cm.  This  resulted  in  reduction  in  the  thrust  produced  to  a  mean  value  of 
-0.32N  and  the  mean  lift  is  0.58  N.  The  minimum  and  maximum  values  of  the  lift  force  during 
the  cycle  are  -0.95N  and  2.3N,  respectively.  Figure  14  shows  the  time  history  of  force 
production  for  this  case. 


Fig.  12.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  20°,  =  1 14 °,/=  3.333Hz. 
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Fig.  13.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  20°,  <t>max  =  114°,/=  1.0Hz, 
rib  spacing  =  2cm  and  with  doubled  rib  lengths. 

As  the  rib  lengths  were  increased  by  a  constant  factor  of  1.5,  the  mean  thrust  improved  to 
-0.25N  and  the  mean  lift  increased  to  1.28N.  The  range  of  lift  force  widened  to  be  between 
-5.68N  and  3.22N.  In  order  to  add  more  curvature  to  the  fin,  the  phasing  of  the  interior  ribs  2-4 
were  modified.  This  resulted  in  an  improvement  in  the  variation  of  the  lift  forces,  but  the  mean 
thrust  reduced  to  -0.46N. 

Another  computational  experiment  was  performed  to  see  the  effect  of  a  more  aerodynamic  root 
section.  For  this  case,  a  NACA0015  section  was  chosen  with  the  maximum  thickness  to  be  1cm. 
The  mean  thrust  remained  nearly  unchanged  from  the  case  shown  in  Fig.  14  and  is  -0.35N;  the 
mean  lift  for  this  case  is  0.95N  with  the  minimum  and  maximum  values  of  lift  at  -0.62N  and 
2.67N.  Therefore,  the  improvement  is  not  significant  to  warrant  a  change  in  the  root  section  that 
will  further  complicate  the  housing  of  the  actuators  for  the  ribs  especially  near  the  trailing  edge. 


Fig.  14.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  20°,  §max  =  114°,/=  1.0Hz, 
rib  spacing  =  1.2cm  and  with  original  rib  lengths. 
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In  order  to  estimate  the  size  of  the  tank  needed  to  test  the  prototype  fin,  an  isolated  fin  was 
immersed  in  a  4'  x  2'  x  4'  computational  tank.  The  results  of  the  simulation  showed  that  the 
force  produced  by  the  fin  was  nearly  unchanged  and  the  walls  of  the  tank  were  not  close  enough 
to  interfere.  This  simulation  did  not  take  into  effect  the  nearness  of  the  free  surface  of  the  water 
which  is  beyond  the  scope  of  this  work. 

The  frequency  of  flapping  was  next  increased  to  2Hz  and  the  force  production  is  shown  in  Fig. 
15.  The  mean  thrust  is  nearly  zero  (-0.06N)  implying  that  the  vehicle  will  be  able  to  sustain  a 
3kt  current.  The  mean  lift  for  this  case  is  0.58N  with  the  minimum  and  maximum  values  of  the 
lift  force  at  -4.55N  and  4.95N,  respectively. 


Fig.  15.  Thrust  and  Lift  production  for  a  UUV  moving  at  3kts,  ar  =  20°,  ())mwi  =  1 14°,/=  2.0Hz, 
rib  spacing  =  1 .2cm  and  with  original  rib  lengths. 

The  pressure  distribution  and  the  force  vector  at  the  instant  when  the  thrust  reaches  a  peak 
during  the  upstroke  and  the  downstroke  is  shown  in  Fig.  16a-d.  During  the  upstroke  high 
pressure  is  observed  in  the  upper  surface  of  the  fin  (Fig.  16b)  and  a  small  region  on  th  elower 
surface  near  the  intersection  of  the  leading  edge  and  the  root  of  the  fin  (Fig.  16a).  Due  to  the 
curvature  of  the  fin  at  this  instant,  a  positive  thrust  is  developed.  During  the  downstroke,  high 
pressure  region  is  present  on  the  lower  surface  of  the  fin,  producing  a  positive  thrust  and  a 
positive  lift  force.  The  coordinate  system  shown  in  this  figure  are  to  orient  the  fin  with  respect  to 
the  inflow  direction  and  is  not  the  actual  location  used  for  computations. 

The  hydrodynamic  power  requirement  for  this  configuration  was  obtained  from  the  force  time 
history  and  the  velocity  of  the  fin.  For  this  purpose,  the  third  cycle  of  oscillation  of  the  fin, 
between  1.35  secs,  and  1.85  secs,  was  chosen.  The  mean  power  required  during  the  cycle  is 
1.573  Watts  which  translates  to  an  energy  consumption  of  0.79  Joules/cycle.  The  inertial  power 
required  to  move  the  fin  can  be  obtained  from  the  moment  of  inertia  of  the  fin,  the  angular 
acceleration  of  the  hinge  rotation  and  the  angular  velocity  of  the  hinge.  The  mean  inertial  power 
required  to  move  the  fin  assuming  unit  specific  gravity  for  the  fin  is  0.1W.  These  power 
requirements  are  shown  in  Fig.  17. 
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Fig.  1 6.  Pressure  distribution  on  the  fin  at  maximum  thrust  production  instants. 


10.000 


5.000  - 


CL 


0.000  - 


-5.000 


■  1 — r  i  i  i  i  i — i  i  i  i  i  ”i  i  ’ 

’  Moan  Power  =  1.573W  —  P„* 

0.600 

“  1 — 1  1  1  1  1  1  I  I  1  I  1  1  1  | 

i  Mean  Power  =  0.0988W  - P  r»mal  ; 

- P.,y 

:  - p-z  : 

■ 

1  1  - P  Total1 

•  1  1  1 

l  i  1  i  : 

i  A '  1 

0.400 

1  1  I  1  J 

;  1  i  . 1  : 

!  W  a.  ! 

| 

■  i  i  A' 

;  i  i  A  V 

:  !  rs\ 

a. 

;  I  '/ill 

i  i  /  i  1  1 

i  i  i 

0.200 

1  II  ||" 

:  \  ill  i  ■ 

:  r  i  i  ‘ 

;  1  l  l  ; 

A 

_ >£ 

~  "l 

<: 

=r*. 

- >-i 

!  1 1  /  downstroke  upstroke' 

!  [upstroke  1  downstroke  1  ! 

•  1  V  \  ysl  1 

i  i  i  i  J  i  i  i  i  i  ■  1  i  i 

0.000 

—1 — 1 _ L.  1  -  1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1_ 

1.25  1.35  1.45  1.55  1.05  1.75  1.85  1.95 

time  (sues.) 


0.00  0.10  0.20  0.30  0.40  0.50  0.60  0.70 

time  (secs.) 


Fig.  17.  Power  requirement  for  moving  the  fin  for  a  UUV  moving  at  V  =  3kts,/=  2.0Hz, 
(a)  Hydrodynamic,  (b)  Inertial. 
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SUMMARY  AND  CONCLUSIONS 

The  3-D  unsteady  computations  of  a  notional  UUV  with  a  flapping  fin  have  been  carried  out. 
Several  parameters  in  the  design  of  the  flapping  fin  have  been  varied.  Based  on  our  previous 
study  on  the  effect  of  rigidity  of  the  Bird-wrasse  fin  on  the  force  production  in  fishes,  the  number 
of  ribs  for  the  current  flapping  fin  is  reduced  to  5.  This  is  done  in  order  to  reduce  the  complexity 
of  the  mechanical  system.  The  important  parameters  that  govern  the  thrust  production  are  the 
angle  of  attack  of  the  root  section  of  the  fin,  the  frequency  and  the  amplitude  of  flapping.  After 
several  parametric  studies  varying  the  angle  of  attack  of  the  root  section  of  the  fin,  it  was  found 
that  20°  was  near  optimum.  The  original  NRL-UM  kinematics  was  also  modified  so  that  the 
amplitude  of  the  oscillation  is  114°and  a  frequency  of  1Hz.  The  rib  spacing  was  set  to  the 
minimum  possible  value  of  1 ,2cm,  from  mechanical  constraints.  This  resulted  in  reduction  in  the 
thrust  produced  to  a  mean  value  of  -0.32N  and  the  mean  lift  is  0.58  N.  The  minimum  and 
maximum  values  of  the  lift  force  during  the  cycle  are  -0.95N  and  2.3N,  respectively.  As  the 
frequency  of  oscillation  is  increased  from  1  Hz  to  2Hz,  the  mean  thrust  is  nearly  zero  (-0.06N) 
implying  that  the  vehicle  will  be  able  to  sustain  a  3kt  current.  The  mean  lift  for  this  case  is 
0.58N  with  the  minimum  and  maximum  values  of  the  lift  force  at  -4.55N  and  4.95N, 
respectively.  The  hydrodynamic  power  requirement  for  this  configuration  was  obtained  from  the 
force  time  history  and  the  velocity  of  the  fin.  The  mean  power  required  during  the  cycle  is  1.573 
Watts  which  translates  to  an  energy  consumption  of  0.79  Joules/cycle. 
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